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R-parity conservation is an ad hoc assumption in the most popular version of 
the supersymmetric standard model. Most studies of models which do allow 
for R-parity violation have been restricted to various limiting scenarios. The 



| single- VEV parametrization used in this paper provides a workable frame- 



work to analyze phenomenology of the most general theory of SUSY without 



R-parity. We perform a comprehensive study of leptonic phenomenology at 
tree-level. Experimental constraints on various processes are studied individ- 
ually and then combined to yield regions of admissible parameter space. In 
particular, we show that large R-parity violating bilinear couplings are not 

ruled out, especially for large tan/3. 
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I. INTRODUCTION 



The Minimal Supersymmetric Standard Model (MSSM) incorporates all of the Stan- 
dard Model (SM) gauge symmetries in its Lagrangian (see for reviews). The MSSM 
Lagrangian is also constrained by a new discrete symmetry, R-parity, defined by 

n = {-\f B+L+2S (i.i) 

where B, L, and S are respectively baryon number, lepton number, and spin. Imposing 
R-parity conservation prohibits baryon number and lepton number violating terms which 
could otherwise lead to superparticle-mediated proton decays on a weak interaction time 
scale, in stark disagreement with observations (for a review and references, see H). Because 
R-parity distinguishes ordinary particles from their supersymmetric partners, the minimal 
model gives rise to a distinctive phenomenology. Supersymmetric particles can be produced 
only in pairs, and the lightest supersymmetric partner cannot decay. These features drive 
many, if not most, supersymmetry (SUSY) search strategies. 

There is, however, no compelling reason to require R-parity conservation. Less restric- 
tive symmetries — conservation of baryon number alone, for example — can be imposed 
to prohibit unwanted proton decay. Furthermore, R-parity is not gauged or required by 
dynamics in the SM or MSSM, and hence there is no theoretical justification for requiring 
its conservation. It is of course possible to devise extensions of the MSSM in which R-parity 
is naturally conserved [|J, but such models remain largely ad hoc. 

Recently there has been a surge of interest in SUSY theories without R-parity. It is 
clear that their phenomenology can differ dramatically from that of the MSSM, and must 
therefore be taken into account in SUSY searches. Introducing R-parity violating terms 
into the superpotentials of these theories complicates any analysis enormously. For that 
reason, most studies make simplifying but otherwise unmotivated assumptions that preclude 
general application of the results. We have adopted a purely phenomenological approach to 
supersymmetric theories without R-parity that provides the framework necessary to greatly 
simplify analyses without necessitating so many a priori assumptions. We introduced our 
approach in 0] and we elaborate on it here. 

The most general renormalizable superpotential without R-parity for a supersymmetric 
model with the minimal particle content may be written as 

W = Sa^QIHtUf + hf^H'Df + hifLIHlEf + ^ H a d H b u ] 

+ e ab[ KjkQiL b jD% + XijkL^LjE k + fi k LlH b ] , (1.2) 

where i, j, and k are family (flavor) indices. The coefficients A and A" are antisymmetric 
in the first two indices as required by SU(2) L and SU(3) C product rules, respectively. The 
SU(2) L indices a and b are shown explicitly contracted with the antisymmetric tensor e a b, 
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with 6i2 = —£2i = —1, to generate SU(2) L singlets. i?-parity conservation corresponds to 
setting Hk = and all A = A' = A" = 0; for baryon-number conservation, only A" = is 
required. 

We can obtain a more compact form for the superpotential W by noting that the SU(2) L 
doublets Lk transform under the same SM gauge group representations as Hd, (SU(3) C , 
SU(2) L , U(1) Y ) = (1, 2, |). (Here the hypercharge Y is normalized such that Q = T 3 + 
Y.) In fact, without the assumption of lepton number conservation, there is nothing to 
distinguish the Hd superfield from the superfields, and they can mix. Their separate 



treatment in Eqn. (L2) is an artifact of starting with the MSSM superpotential and adding 
R-parity violating terms. Therefore a more appropriate form for W, with extended flavor 
indices a and (3 running from to 3, is 

W = e ab [h^H b u Uf + \' iaj Q?L b a Df + X af}j L a a L b p Ef + pL a L a a H b u ] + X^D-f 'f)fU° . (1.3) 

£7(4) flavor rotations can be used to transform between bases of the four superfield doublets. 

It is in principle possible for the neutral scalar component in each of these doublets to 
acquire a non-zero vacuum expectation value (VEV). The key to our approach, the single- 
VEV parametrization J5],|6|], is to rotate the doublets into a basis in which Lq alone bears a 
non-zero VEV. The remaining admissible leptonic flavor rotations are sufficient to diagonal- 
ize the lepton Yukawas h\ k = 2X i0 k = —2X ik which are then given by ^diag{m 1; m 2 , m 3 }. 
There is then no additional freedom to set the fi^ bilinear coefficients equal to zero; to 
maintain complete generality they must be left arbitary. 

As implied above and discussed in care must be taken to specify what choice of 

flavor basis (if any) is implied when a specific set of of RPV parameters is given. For 
example, if the sneutrino VEV's and ///- bilinear terms are left arbitrary, then they are not 
truly physically independent because of the freedom to rotate between bases. In another 
extreme, setting all sneutrino VEV's and bilinears to zero results in a loss of generality, 
while still not being sufficient to uniquely determine the flavor basis. The common approach 
of using i?-parity conserving MSSM particle states with the ad hoc addition of a few RPV 
trilinear couplings results in ill-defined RPV parameters (since such parameters are in general 
basis-dependent and the flavor basis is not specified). In principle, that renders the analysis 
internally inconsistent. In practice, this approach can reasonably approximate some regions 
of parameter space, but an inherent ambiguity remains. One way out of this quandary 
would be to construct "basis-independent" observables. This has been explored with strong 
a priori restrictions placed on the sneutrino VEV's || . However, the observables were found 
to be phenomenologically "messy" and impractical for most experimental situations even in 
this limiting case. The alternative resolution is to carefully choose a convenient basis that 
renders the experimental analysis as simple as possible — the path taken in this work with 
the single- VEV parametrization. 

As we described in 0, this parametrization has the advantage that the tree-level mass 
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matrices of all fermions in the theory are independent of all trilinear RPV couplings. In par- 
ticular, R-parity-violating contributions to leptonic phenomenology at tree level are almost 
entirely determined by the ///.. In this paper, we examine this tree- level leptonic phenomenol- 
ogy in detail within the single VEV parametrization, focusing on the constraints that can 
be obtained from existing experimental results. The rich variety of relevant data include 
precision measurements at the Z° pole, charged lepton and pion decay rates and branching 
ratios, neutrino mass bounds, lepton-neutrino scattering cross sections, and limits on neu- 
trinoless double beta decay. A complete list of the measurements we will use can be found 
in Table I. 

In the present analysis, we assume that these processes are mediated by gauge bosons, 
i.e., on-shell or off-shell Z°'s or W ±J s. Scalar intermediaries are also possible for various 
processes, and these would re-introduce dependence on the trilinear RPV couplings (and 
other RPV parameters from the scalar sector) even at tree-level. We assume that contribu- 
tions from such scalar modes are negligible due to kinematic supression from the presumably 
much heavier (relative to the gauge bosons) scalars. It is certainly possible in principle that 
relative strengths of the couplings involved and not-quite-so-heavy scalars may conspire to 
invalidate our assumption. However, for these to affect our experimental bounds (or for 
that matter, to be larger than the loop corrections we omit), one has to invoke very special 
and/or unlikely combinations of parameters. Including the additional parameters needed 
to accommodate these tiny corners of parameter space needlessly complicates the analysis 
and will be avoided here. In what follows we will point out neglected scalar modes where 
appropriate. 

We begin our analysis in Section II, where we examine the neutral current interactions 
of the charged fermions (charged leptons and charginos); in Section III we do the same 
for the neutral fermions. We then consider charged current interactions in Section IV. 
Having considered the constraints separately, we then combine them in section V, where we 
determine exclusion regions in the parameter space from both charged and neutral current 
processes. Finally, we present a summary of our results and our conclusions in Section VI. 

II. COLOR-SINGLET CHARGED FERMIONS 

The color-singlet charged fermions, the charged leptons and the charginos, are linear 
combinations of two-component Weyl spinor charged winos (written as \±), the spinors of 
the Ei superfields, and the spinors for the charged superfields of the L a and H u SU(2) L 
super doublets. In terms of Dirac spinors for a theory with spontaneous symmetry breaking, 
denoting spinors corresponding to the superfields in Eqn. (|1.3p by the corresponding lower- 
case letters (and where A is the wino and ip is the "higgsino" built from Iq and h u ), the 
contributing Lagrangian mass terms are given by 
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Here the two vacuum expectation values are in general complex, as are the /x Q 's and the 
soft gaugino mass M 2 . As is customary, all of the parameters entering the fermion mass 
matrices will henceforth be assumed to be real. Potentially very interesting CP- violating 
effects related to the possible complex nature of these parameters will be relegated to a future 
paper. On the other hand, it must be strongly emphasized that there is no loss of generality 
concerning the trilinear RPV terms — these need not be set to zero. They simply do not 
appear at tree level in the fermionic mass matrices in the single- VEV parametrization. 

We now wish to find the mass eigenstates. In the MSSM, the Yukawa mass entries, 
mi, m2, and 777.3 i m i — hfi v d in the single- VEV basis), are identical to the known physical 
charged lepton masses. Without R-parity conservation this correspondence is spoiled by 
the presence of the fa in M.q- Masses of the well-known charged leptons (and in fact the 
eigenvectors and mass eigenvalues for all five physical states) now result from diagonalizing 
the 5x5 chiral mass matrix M. c : 



UlM c U R = diag{M cl , M c2 , m 1 , m 2 , m 3 } 



(2.3) 



where Ul and Ur are unitary matrices which diagonalize M. c , and m e = mi, = m 2 , and 
m T = 7773. In general, the rrii depend on the and the fa. 

Therefore for non-zero fa the input Yukawa parameters for M.q need to be determined, 
i. e. for a fixed set of fa, we need to find the values of mj that give the correct physical masses 
rrii. This is done by writing a system of linear differential equations for the infinitesimal 
change in the m^'s due to an infinitesimal change in the /ij's. Beginning with all fa = 
(where the solution of the system of equations is known), an acceptable solution for a chosen 
set of fa's is then obtainable via numerical integration HlO| . Note that only the three 'lepton' 
masses need to be fixed in this way. The heavier so-called 'chargino' masses M cl and M c2 
also depend on the fa, but we do not (yet) have experimental constraints telling us what 
the physical chargino masses should be. 

We now consider interactions with the Z° boson, following Ref. but using the single- 
VEV parametrization. The couplings of the five mass eigenstates to the Z° boson are given 
by 
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X+X-Z° 
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where xt 



A L 



2 cos 6* u 

Xa = (XuX2j e JixJr), Pr.l = |(1 ± 7s), and 

5 ab (l - 2sm 2 9 v 



Ul a Ul b + 5 ab (l - 2sinX) = 5A L ab 



A R ab = 2£#£/# + U 2 R a U 2 R b + 26 ab sm 2 9 w = 5A R ab + 25 ab sin% 



(2.4) 



(2.5) 
(2.6) 



(notation here follows that of [TT[] except for the ordering of basis vectors). The 5 ab terms 



ab 



are the SM expressions, and deviations from SM expectations originate from non-zero SA 
and SA R b . The anomalous coupling of any two charged fermions to the Z° can thus be 
determined in terms of the Ul and Ur matrices found numerically from diagonalizing A4 C , 
and this is precisely what we do to obtain the parameter space plots presented later in this 
paper. The exact analytic expressions for the eigenvalues and eigenvectors of the 5x5 Ai c 
prove cumbersome and not very illuminating and we do not reproduce them here. However, 
to improve our understanding of the results of the exact numerical analysis, it is very useful 
to consider a couple of interesting analytic approximations. 

First, we treat the R-parity violation as a perturbation, taking the limit in which the ^ 
are small. If the /ij were zero, the (now MSSM) chargino sector, which is the upper 2x2 
portion of Ai c , would be orthogonal to the already diagonal SM charged lepton sector. So 
first we introduce new 2x2 rotation matrices R L (9 L ) and R R {6 R ) such that 



/ 



M 9 



V 7^ 9Vd /io 



R R = diag{M cl ,M c2 } 



(2.7) 



M cl and M c2 are the chargino masses in the /ij 
M. c we find 



/ 



M c ' = RlM c Rr 



Mr 



MSSM limit. Applying this rotation to 











M c2 

/ij sin^R /ij cos9 R m 1 

fi 2 sin9 R /i 2 cos9 R rri2 

H 3 sin9 R fi 3 cos9 R m 3 



(2.8) 



To obtain the general 5x5 Ul and Ur matrices from the MSSM's 2 x 2 R L and R R matrices, 
we treat the 3x2 off-diagonal block containing the perturbation. Then 



Ul,R 



P\,r Pl,rVl,R 



L.R 



'3x3 



(2.9) 
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where the elements of the 3 X 2 Vl and Vr matrices are given by 



T/il _ /o fijM w cos/3 -sril _ fnmi HiM w (M2 sin/3+/2 cos/3) 

V L - V R - V^ Mq Mq 

T/ i2 _ _\hM 2l T /i2 _ m, lh(M$+2M% v cos 2 I3) 

V L - ~Mj~ V R - M M% 



(2.10) 



and where 



Ml = n M 2 — sin2/3 (2.11) 



which is in fact the determinant of the 2x2 MSSM section of Ai c - In the limit of large 
tan/3 and/or /i M 2 3> M^,, M 2 tends toward fJ, M 2 . 

In terms of the Vl and Vr matrix elements, we can now write the deviations of the 
fermions' coupling coefficients to the Z° boson from the SM case as 

ifc=V?Vf = 2^1^l (2.12) 



ab 

and 



Ml Ml 



5A R ab = 2V R 1 Vi 1 + VgV* 
_ ///,///, Hiiij 



Ml Ml M 4 



[M 4 + 4M 2 .(M 2 + /i M 2 sin2/3 + /i 2 cos 2 /? + il^cos 4 /?)] , (2.13) 



for {a, i}, {b,j} = {e, 1}, {/i, 2}, {r, 3}.|] This provides us with simple, quantitative expres- 
sions of how the gaugino and higgsino contents of the charged lepton mass eigenstates (i.e., 
of the e, the /i, and the r) affect the Z° leptonic decay widths. These effects, present when 
the Hi are non-zero, can lead to non-universality among the Z° leptonic branching ratios or 
to mixed-flavor Z° leptonic decays. 

Notice that the deviations of the Z° coupling coefficients are in fact proportional to ^jt, 

M o 

which must be small in this perturbative approximation (this tells us quantitatively what 
it means for the /i« to be "small"). In addition, the 5A% b are suppressed by the factor m ^ J . 
Thus one might expect that the effects of the b~A L ah will dominate those of the 5A^ b . However, 
the 5A^ b are themselves proportional to cos 2 /? which will strongly suppress their values in the 



In the formulae above and those that follow, index labels a and b are reserved for the physical mass 
eigenstates, a,b = e, fj,,r (or the heavier charginos). Index labels i and j are used for basis state 
parameters in the single- VEV parametrization, i,j = 1,2, or 3. In the small-fii approximation, a 
and i (or b and j) have a simple one-to-one correspondence. Hence, equations could be written 
with just one pair of index labels. However, both pairs will be kept to make clear that this is not 
true in general. 
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large tan/3 regime (for example, cos 2 /? ~ 10 3 for tan/3 = 45). Neglecting terms proportional 
to cos/3 in the 5A^ b (recall sin2/3 is small when cos/3 is small), we find 



l (^ i ,^W^ M 2 [(M 2 /M w ) 2 + 4] 



1 



tar?/3) 



M n 2 



-^-[(M 2 /M w ) 2 + 4] (large tan/3 or /i M 2 > M 



(2.14) 



Thus if M 2 is large and |/x | is small, the effect of deviations from the SM for the right-hand 
component of the Z°-t-t coupling (and even the Z°-t-]i coupling) can be as significant as 
those for the left-hand component. 

The Z° branching fraction into any pair of colorless charged fermions (with M z > M Xa + 

M x b ) is g iven b y 



Br(Z 



o 



a 2 



M| ' Af| J M 7 , 



24cos 2 #„ 



M la+ M l b ( M la- M lj 2 , 6M Xa M Xh Kb A ab 



Alb - 



.4 



ab 



where a 2 = g\j^K and 



(2.15) 



12 _ I |2 , I AB |2 



ab I 



(2.16) 



The experimentally determined total decay width of the Z° is V z = (2.4948 ± 0.0075) GeV 
||12|| . The kinematic A-function is A(a, b, c) = (a — b — c) 2 — Abe. In the small fa perturbative 
approximation, and for tan/3 not too large, we have for the leptonic Z° decays, 



A\ ~ (1 - 2sinX) 2 + (4sinX) + 2(1 - 2sin 2 ^)&4 J I 

Ml cos 2 /3 / iu ^ 2 
= .5027 + 2.148 w - f — 
Ml \M 



and 



.4 



ab 



4M 4 cos 4 /3 (jh_\ 2 



(o,i ^ 6, j) . 



(2.17) 



(2.18) 



If tan/3 is large, then the right-hand component cannot be neglected, as noted above, and 
the A 2 and A 2 formulae should be modified to 



and 



At 



A 



1 - 2sinX) 2 + (4sinX) + 2(1 - 2sm 2 9 w )8A L TT - 4sinX&4* T 



.5027 + 



\K T \ 2 



211g M 2 cos 2 /3 _ >926 ^|(M 2 + 4M|)M : 



M 2 





4M 4 cos 4 /3 + m 2 m 2 (M 2 + 4M 2 .) 2 M 



M 4 



M 12 



Ma \ 2 / N ~ 



MJ \M 



(2.19) 



^2.20) 
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Note here that the numerical value sin 2 6L 



0.2315 fl2H has been used. This effective value 
couplings. Note also that the 



for sin 2 6 w absorbs the SM radiative corrections to the Z°£, 
two terms inside the bracket in Eqn. (|2.19| ) enter with opposite signs, implying that in some 



region of the parameter space, the deviation of T(Z C 



7 7 



from the SM prediction could 



be suppressed by cancellation between these two terms even if [i 3 is quite substantial. 

We now apply these results for couplings and branching ratios to specific processes in 
order to obtain limits on the [ii. 



A. Mixed-flavor leptonic Z° decays 

For non-zero /Vs, the Lagrangian of Eqn. ( |2.4| ) leads to the tree-level flavor-violating Z° 
decays, Z° — > efi, er, and fir. The predicted branching ratios for these decay modes are 
given by Eqn. fl2.15|) . The experimental bounds from LEP on these processes are shown 



in the third column of Table I. In the small-/Zj approximation, Eqn. ( [2.1 8|) translates these 
constraints into the following bounds: 

^<Ml + tan^ (2.21) 

lvl o lvl w 



£12 


= 1.8 x 10~ 3 


from Z Q - 


4eV 


£13 


= 4.3 x 1(T 3 


from Z° - 


-> e ± r T 


£23 


= 4.7 x 1(T 3 


from Z° - 





Note that bounds apply to the products \[ii[J.j\- That the constraints can be cast in such a 
simple form and in terms of so few RPV input parameters is a key strength of the single- VEV 
parametrization. Note also that the constraints become weaker as tan/3 increases. 



B. Flavor- violating charged lepton decays 



The Z £i£j couplings can also produce the tree-level FCNC decays of [i and r via a 
virtual Z°. These branching ratios are given by 

BHt- - ItKK) = (^)* ^ {%Jl + I 2 + , (2.23) 

where Ti a is the total decay width and m^ a the mass of the decaying £ a (masses of the 
daughter leptons are neglected). It is assumed the virtual intermediate is an off-shell Z°; 
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possible slepton/Higgs intermediates are assumed to be heavy and therefore to yield negli- 
gible contributions. Note from Eqn. fl2.23| ) that when the odd-flavored daughter lepton 
has the charge opposite to the parent lepton, the branching ratio is suppressed by an extra 
factor of A 2 ab (a ^ b) compared to the case where the odd-flavor daughter lepton has the 
same sign as the parent. Actually, the latter case also has a contribution analogous to Eqn. 
( |2.23| ); however, this is insignificant relative to the Eqn. ( |2.22j ) contribution. 

The experimental limits on these decays are again given in the third column of Table I. 
Using Eqn. ( [2.18|) along with the lowest order value for A 2 CC (= 0.5027) from Eqn. ( p.l7|) , 
bounds can be cast in the same form as in Eqn. ( 2.21| ), with : 

/C12 = 1.4 x 10" 6 from \T e~e + e~ 
/C13 = 4.5 x 10~ 3 from r~ — > e~ ^ \T 
JC 12 = 4.3 x 10 -3 from r~ — > iTe + e~ . 

We see that while the r decay constraints are comparable to the Z° — > t£ [t = eor/i) 
constraints, the fi — *> eee constraint is much more strict than that for Z° — > efi due to the 
much stronger experimental bound on fi — > eee. In Fig. 1, predicted values from the exact 
numerical calculation of this branching ratio are shown as a function of the \ii for the MSSM 
parameters M 2 = fi = 200 GeV and tan/? = 2, 45. The constraint is very stringent for small 
tan/? and remains relevant, although much weaker, even for large tan/5 (unlike the other 
constraints discussed so far). Plots for Br(Z° — > e\i) would look very similar except with a 
much weaker experimental bound. Note that in Fig. 1, the useful RPV parameter 



M« = yjtf + tf + til (2.24) 

is introduced to permit different [i± : /i 2 ratios to be plotted simultaneously. This parameter 
will appear repeatedly in subsequent discussions. Interpreted as a constraint on fi 5 , Fig. 1 
shows that the fi —> eee constraint can be evaded by supposing a strong hierarchy among 
the /ij's (i.e., fi r <C fi 2 /i 3 ), as can also be clearly seen from the approximate Eqn. ( |2.21| ). 

Fig. 2 shows Br(r~ — > e~fi + fi~) for tan/3 = 2, again based on exact numerical calcu- 
lations. There is no meaningful constraint for tan/3 = 45 due to the much weaker experi- 
mental bound. We need not explicitly show results for the remaining processes: results for 
Br(r~ — > e~e + e~) and Br(r~ — > ^ fi + ^) are almost identical to Fig. 2, as are those for 
Br(r~ — > /U~e + e~) and Br(r~ — > if the roles of fi 1 and fi 2 are interchanged. Plots 

for Br(Z° — > er) and Br(Z° — > fir) also yield very similar results. Note that in using Eqn. 
( |2.21| ) to obtain constraints on r decays in the small /z; approximation, extra contributions 
from A^ T have been neglected. These A^ r contributions are only significant at large tan/3, 
where constraints on \ii from this process as given by Eqn. ( 2.21 ) are superseded by limits 
from other processes. This is confirmed by exact numerical results. Finally, note also that 
for r decays in which the odd-flavored daughter lepton has the flipped charge, the result is 
proportional to A ejl and thus very strongly limited by constraints from Br{t~ — > i^^t^c) 
such that Br(r~ -> £+£~£~ )<10~ 18 . 



10 



C. Universality violations at the Z° peak 



Eqns. ( p. 17] ) and ( |2.19| ) also produce deviations from SM predictions for Br(Z° 



which can break lepton universality. Fig. 3 shows leptonic partial widths as a function of 
/i 5 for several choices of : fi 2 : /i 3 , again using the representative MSSM parameter point 
M 2 = /i = 200 GeV and tan/3 = 2. The experimental 3a bounds are only exceeded for 
fi 5 values in excess of 50 GeV. As the value of /x 5 increases, the small-/ij approximation 
loses its validity. The numerical results show that for sufficiently large //j values, the partial 
decay widths stop increasing as the /ij's increase and in fact turn over and decrease. This 
behavior is common to all the decay widths discussed so far (although the maxima will in 
general occur at different ^5 values, sometimes above the upper limits on the plots shown). 
This behavior can be understood via the "large //j" approximation employed in subsection 
E below to derive chargino masses. Note that if, at some point after reaching its maximum, 
the deviation in a partial width again drops below its experimental limit, then the large /I5 
values above this point are again acceptable. Thus it is possible that no upper bounds can 
be placed on /x 5 from these processes, but instead either only a finite range of /is values is 
excluded or (if the maximum is too low) no /i 5 values are excluded at all! For further details, 
see fC3"| . In practice, other constraints will rule out arbitrarily large /i 5 values. 

Universality constraints on Z° — > £ a £b decays are quantified via the observable uj£ 



MiM = ^^CC)-r(Z^^1 = Al - A\ h Qj - ^)Mj cos^/3 



where the first equality follows from Eqn. ( |2.4| ) and the second from Eqn. (|2.17|) . LEP 



experimental measurements on the partial widths can now be translated into the restrictions 
on the U^. atb ^ listed in the third column of Table I. These are all compatible with the SM 
prediction of U^ aib ^ = 0, and, neglecting the (small) nonzero central values, translate into 
the following bounds on the fif 

1 — M 2 

— ^-^^^•(l + tan 2 /?)^ . (2.26) 

The XVs have values of £ 12 = 2.05 x 10~ 3 , K 13 = 2.33 x 10~ 3 , and £ 23 = 2.62 x 10" 3 ; they 
are comparable to the /Cy's of Eqn. ( |2.21 ) obtained from flavor- violating Z°-decays and r- 



decays. For large tan/3, deviations from the SM are highly suppressed; i.e., very high values 
of // 5 (well beyond the range of validity of this approximation) are allowed. Exact numerical 
studies also confirm that these constraints vanish. Using the above formula from the small 
/ii approximation, we find that yU,'s (or more precisely their difference in magnitudes) as 
large as M become allowable for tan/3 ~ 20. In fact, for large tan/3, Eqn. ( |2.19| ) should be 
used for the Z° — > t + t~ partial width. Cancellation among terms in this equation would 
further weaken any surviving Uf r ar) bound on /i 3 . 
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D. Leptonic left right asymmetry 



Predictions for left-right asymmetries in Z° leptonic decays, which are defined by 



A 



\A L I 2 



\A R I 2 
I Awl 



I A L I 2 
I Am I 



+ \A 



R 12 
aa I 



(2.27) 



follow immediately from Eqns. ([2. 171) an d ( [2.19|) : 



A = 4 SM) + 4.273 M - COS ^^^ V " 



Ml 



_) [a = e , j = 1 or a = = 2) 



(2.28) 



Mjcos 2 /? m 2 (M 2 + 4ilQlf 2 

~" > ^ """" 



M 



^2.29) 



First note that the |/ij|'s enter individually rather than in products or differences of two 
distinct |/x;|'s. This makes the left-right asymmetries potentially very useful in distinguishing 
effects from the three /ij's. Note also that non-zero /Vs always increase the AJs from their 
SM values, and that the SAf T contribution now reinforces that of 5A^ T . This could be 
important if these contributions cancel in yl 2 r which enters into the previously-discussed 
RPV r effects. Another immediate consequence is that for the favored case of /i 3 > fi 1 , 

Eqns. ( 2.28 ) and Q2.29 ) and conclusions drawn from them are valid in the small-/ij ap- 



proximation. In Fig. 4, the asymmetries are shown using exact numerical calculations (again 
for M 2 = /i = 200 GeV). For large values of // 5 , the ^4 a 's cease rising with increasing /x 5 , 
deviating from the approximate behavior of Eqns. fl2.28|) and (|2.29|) . As with the Z° leptonic 



partial decay widths, a maximum is reached and then the slope turns negative. For low 
tan/?, this occurs for /i 5 values excluded for other reasons. But for high tan/3 it may occur 
in admissible regions of the parameter space. Also notice from Fig. 4 for tan/3 = 45 that A T 
shows an increase relative to A e and A^ arising from the right-handed contribution. 

Eqns. ( |2.28j ) and (|2.29| ) require A[ S as input — this must be determined independent of 



the leptonic asymmetry measurements. Unfortunately, Af depends strongly on sin 2 ^; for 
instance, Af M) decreases by more than 11% when sin 2 ^ is increased by 1%. The effective 
value of sin 2 ^ for the Z°££ coupling depends on radiative corrections |15|]; sin 2 ^ = 0.2315 
employed here yields A ( f M) = 0.147, but includes only SM corrections and not additional 
corrections depending upon SUSY parameters. Thus uncertainty in the effective value of 
sin 2 #„, leads to even larger uncertainty in A\ , which also includes beyond-SM contributions 
not apparent in the simple separations seen in the formulae above. In an attempt to reduce 
such uncertainties, Ref. [|TTJ (following the idea put forward in Ref. |16[) suggests using 

AA* , ^ (2.30) 
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instead of the individual AaS. AA a b = in the SM. In the small fi{ approximation, where 
A a — Af M) is also small, 

AA b = - l) Ul aih) (2.31) 

(as noted in Ref. [jlll), if 5A^ a 's are neglected. If this is permissible, then the AA a b con- 
straints have the same dependence on the //j's as the constraints for the U^. a£b \ providing 
a compatibility condition for the RPV framework. The M£ a 's may safely be neglected for 
the first two generations, meaning the allowed window for AA efl given in Table I follows 
immediately from that for (though perhaps a small amount of extra widening to re- 

flect the uncertainly in A^ M) should be included). Unfortunately, the present experimental 
situation is far too imprecise for any incompatibility to be seen. As noted previously, 8A^ T 
may not be negligible, especially for larger tan/3. Adding these terms to AA mT and U b ^ mT ' > 
means that the value of one is no longer determinable given only the value of the other. Fur- 
thermore, deviations from Eqn. ( |2.31|) , which predicts |A^4 mT |<0.03, indicate a substantial 
contribution from 5A* T (within the framework of the small-/^ approximation). 

In the t case, then, since Uf£ experimental constraints cannot place bounds on the 
A^4 ar 's, direct At measurements must be considered. These are listed in Table II and used 
to derive the 3a bounds on the A.4. aT 's shown in Table I. The A e and A T used in the AA eT 
bounds of Table I are obtained from LEP measurements of r polarization (note however that 
the SLD group at SLAC (T7J measured A e directly with their polarized electron beam and 
obtained a substantially higher value for A e ). The AA^ T bounds are obtained from LEP 
forward-backward asymmetry measurements (which provide the best A^ value but larger 
uncertainty in At). Note that the bounds are considerably weaker than the 0.03 value 
given above, meaning that substantial deviation of A T from Af M) including an important 
contribution from SAf T is possible. 



E. Limits on charginos 

We now turn to constraints associated with the remaining color-singlet charged fermions, 
the charginos. The term "chargino" is here applied to the two (heaviest) mass eigenstates 
remaining after the other three eigenstates in Eqn. ( |2.3| ) are fixed to give the "leptons" with 
well-known experimentally observed masses (mi = m e , m 2 = m M , m 3 = m T ). In the MSSM, 
the charginos have masses given by M cl and M c2 in Eqn. (|2.7| ). For nonzero values of the /ij, 
the chargino masses are modified to M cl and M c2 in Eqn. Q2.3|) . If the \L{ are "small," then 
M a — M cl and M c2 ~ M c2 ; however, larger values of the /i^s can produce more pronounced 
effects. From Eqn. (|2lf), 
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Mc(M c 



gvu. 
gvu.. 



Now by applying the rotation 
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where i? 
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\0J 



[2.32) 



and /x 5 is given by Eqn. fl2.24| ), M. C (M. C )^ can be rotated into the form 



£A + M 2 2 2| M 2 + ^ ^/i 5 
3MM 2 + ^fi ^+/i 2 w s 



V2> 



gvu. 


















+ 







2x2 



Ri 



m\ 








\ 



m\ 
V ml) 



Re 



(2.33) 



The matrix R 5 will appear again in the next section when the color-singlet neutral fermions 
are considered. In the limit where // s 3> m 3 ~ m T (referred to henceforth as the large 
/ii approximation^) , the second matrix of expression (|2.33|) that is proportional to the m/s 
may be dropped. This leads to simple analytic formula? for the chargino mass eigenvalues: 



m: 



Ml + 2M 2 W + /i 2 + n 
1 



± 



2 ( a i + "2) ± 2 



(/x 2 + /i 2 - M 2 - 2M 2 cos2/5) + 8M 2 (M 2 sin/5 + /i cos/3) 2 
1 



1/2 



a 2 ) 2 + 2(gv u M 2 + gv d n ) 



M 2 C1 + M 2 + £ 



± - 
2 



(M, 



M 2 -/i 2 )~-4/i 2 cos 2 ^(M ( 



M 2 



1/2 



(2.34) 
(2.35) 
(2.36) 



where 



Note that this is not the converse of the small /Xj approximation, which requires /ij <C Mq. 
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a- = ^ + yU 2 + /i 2 an d a 2 = ^ + M 2 . (2.37) 

Within this approximation the chargino masses only depend on the fa through the single 
parameter fa. The roles of fa and fa in the mass formulae are very similar. Note from Eqn. 
( p.36| ) that the chargino masses reduce to the MSSM chargino masses when /i 5 = 0. Also 
note that a non-zero value for fi 5 increases the lighter chargino mass; hence, some region of 



the M 2 - fa parameter space that is ruled out in the MSSM by the chargino mass bound (18 
can be re-instated when fa ^ 0. However, correct chargino mass bounds for RPV scenarios 
require analysis of the decay modes. One such analysis was recently performed for a specific 
bilinear RPV model fT9| |; a general analysis in the single- VEV parametrization, including 



contributions from both bilinear and trilinear RPV couplings, is currently in progress. 

It is possible to simultaneously have /i 5 3> m T and fa <C M , meaning that both the 
"large-/ij" and the "smaH-/ij" approximations will yield quite accurate results. In fact, results 
from the two methods agree well with each other and with the exact numerical results in 
most relevant regions of the parameter space. Using the exact numerical calculation, the 
region in MSSM parameter space that would be excluded by a 90 GeV chargino mass bound 
is plotted in Fig. 5, from which it is clear that constraints based on such an analysis will be 
significant if M 2 and/or fi are small. 

Bounds from Z° decays into two charginos or a chargino and a lepton follow from the 
general formula for Z° decays into a pair of charged fermions, Eqn. ( |2.4| ). However, de- 
termining the exact values for these bounds requires detailed information about chargino 
decays from detector simulations which is beyond the scope of this work. Therefore, only 
a conservative bound is applied: Z° decays involving at least one (on-shell) chargino are 
required to have branching ratios of less than 10~ 5 (as noted in Table I). 

Turning to the lepton sector, we see that in the limit where mf = 0, the matrix J\A C (A4 C Y 
of Eqn. (|2.32| ) [or Eqn. ( [2.33| )1 has three zero eigenvalues. Thus, in this approximation, the 
masses of the physical leptons are zero. Despite this degeneracy, the correct eigenvectors can 
be found by requiring that they do not differ significantly from the exact massive eigenstates 
obtained when the non-zero m ; are retained if /x 2 is large. Three such massless eigenvectors 
for A4 c (Ai c Y with m\ = are given by: 



^2.38) 



(the Aj ee yj M ( 4 + a 2 /i 2 are normalization constants). Since these are degenerate states, 
any linear combination of the three will be a massless eigenstate. The natural hierarchy of 
the rrii values (m 3 > m 2 > mi) requires the following choice of orthogonal eigenvectors: 
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-M 2 fa, 


Ml 
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-M 2 fi 2 , 
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Ml 


f 


/A 3 



15 



|4 > oc \i\ > 

= , -M^ x Ml , , f /A e 

= (^/i 2 , -M 2/ i 2 , -a^^/M* , A 2 /M 2 , f M 2 /(A e A /i ) 

% > ex |f 3 > -\£[ >< e[\£' 3 > -\e 2 >< £' 2 \e' 3 > 

= (fif// 3 , -M 2f i 3 , -a^/ijM* , -a 2 /i 3 /i 2 /M 2 , A 2 /M 2 ) T M 2 /(A M A r ) 
where the normalization constants are: 



(2.39) 



A e = yX 4 + a 2 tf , A M ee ^M 4 + a 2 (/i 2 + ,u 2 ) , A r ee ^/M 4 + a 2 (/i 2 + /i 2 + /x 2 ) . (2.40) 

The first components of these three vectors are the elements Uf} of the left rotation 
matrix in Eqn. (|3|) : 

ut = 9^}H vpmpJfi mi = ^/!X . (2 . 41) 
L v/2 A e L v/2 A e A M v/2 A M A T v ; 

These quantities appear in the expression for the anomalous coupling coefficients SA^ b of Eqn. 
( p.5| ). Note that in the limit /ij — > 0, A a — > Mq and we recover the "small /V approximation 
results given by Eqns. ( |2.1(J| ). The dependence of the normalization constants (the A a 's) on 
the /ij's indicates that the decay widths will decrease (for /i and r) as the /ij's increase, 
and exposes the reason behind the maxima seen in Figs. 3 and 4. Thus if an experimental 
bound actually allows the maximum value (for instance at large tan/3 where there is strong 
suppression), then that bound may provide no constraint on fi 5 (or equivalently on the ytij's) 
at all. Even if the maximum is excluded, a window of large (possibly very large) /ij's may 
be allowed. However, this region of the RPV parameter space is largely ruled out by limits 
on neutrino masses and bounds from neutrino scattering and charged current processes. 
Therefore we do not consider very large /ij values further at present and instead proceed to 
a discussion of these additional constraints. 



III. COLOR-SINGLET NEUTRAL FERMIONS 



In the single- VEV parametrization, the Lagrangian terms contributing to the color- 
singlet neutral fermion (neutrino and neutralino) masses may be written as 

C-3 -L\C t\c h u oc h d oc v c v c v G^M M ^-i\ Q -i\ 3 h u ° tif v x u 2 z/ 3 ) + c.c. (3.1) 

where h u ° , hf , and Vi (i = 1-3) are the Dirac spinors associated with the neutral superfields 
in H u , L , and Lj, respectively, of Eqn. ( |1.3|) and — i\ and —1X3 are respectively the bino and 
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wino components. C refers to charge conjugation acting on each spinor. Since neutrinos are 
assumed to have zero Yukawa masses, the same -R5 rotation that was used when the Yukawa 
masses for the charged leptons were neglected is also applicable to M.^. Now there is no 



approximation, no Yukawa part to throw away; two massless neutrino states decouple |20 
leaving a 5 x 5 matrix: 
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(3.2) 



A. Neutrino mass (tree-level) 

The single massive neutrino that results from diagonalizing this tree-level mixing matrix 
was discussed in the first paper in this series Here again we denote the massive state by 
1^) — Wi) + 1^) + ^ l 1 ^) ■ Approximate analytic formulae for its mass were found to 
be: 

1 tfv 2 cos 2 (3 (xg 2 + g' 2 ) 

m " 5 2 /i [2xM 2 fi - (xg 2 + g' 2 ) v 2 sin/3 cos/3] 1 ' ' 

(where v 2 = v 2 + v\ and M x = xM 2 ; x — | tan 2 ^ assuming gaugino unification, as is done in 
all numerical calculations) from a "seesaw" approximation in which /z 5 is taken to be small 
and 

1 /j, 2 v 2 cos 2 {3 (xg 2 + g 12 ) 

from a perturbative treatment in which the electroweak symmetry breaking terms are re- 
garded as small but the magnitude of /i 5 is not restricted. Key features to note are: 

• The mass has a simple dependence on only one RPV parameter, fi 5 . This is to be 



contrasted with results found without using the single- VEV parametrization The 
formulae above are also quite general. In particular, the trilinear terms have not been 
set to zero — they simply do not contribute to the tree-level mass formulae. They will 
reappear at the one-loop level; however, loop effects are expected to be small compared 
to those at tree-level [p2j . 
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The approximate expressions, Eqns. ( |3.3|) and fl3.4j) , indicate that the neutrino mass 
is proportional to n\ if fi 5 <C fi . Eqn. Q3.4Q also shows that for /x 5 3> /x , the neutrino 
mass approaches a constant asymptotic value. The features of these approximate 
formulae are confirmed by exact numerical calculations, as shown in Fig. 6 for the 
MSSM parameters M 2 = fi = 200 GeV and tan/3 = 2, 45. Note in particular the linear 
rise in logm^ seen in the log-log inserts for low log/z 5 values, as expected if m U5 oc p? h . 
Note also that there is no decrease in m U5 for high /i 5 values, in contrast to the case 
of several constraints in the previous section. Therefore, unless the asymptotic value 
falls below the experimental bound, the \x h upper limit limit from m us will close any 
large-/i 5 window. 

There is strong suppression both at high tan/3 due to the cos 2 /3 factor^ and at high 
M 2 . Inverting Eqn. (|3.4|) yields a bound on /i 5 : 



;/2 < 4x/i M 2 m„ M 

^ <; W/3 (xg* + g'i) - AxM 2 m UM ' 1 ' } 

clearly demonstrating the high tan/3 suppression. The exact numerical results in both 
Figs. 6a and 6b and Figs. 7b and 7d also illustrate that the asymptotic limit is much 
lower for large tan/?. In addition, as M 2 increases, the denominator goes to zero, 
beyond which /x 5 is unconstrained. 

If we assume // 5 = // 3 4=> (fi 1 : /i 2 : ^3) = (0:0:1), then the recently-improved 
mass bound on v T from LEP, v T < 18.2 MeV (shown by solid horizonal lines 
in Figs. 6ab), implies that, for tan/3 = 45, there is no m Vr constraint on /i 5 when 
M 2 >280GeV (Fig. 7b). Analogous arguments can be made for /x 5 = /x x and \i h = // 2 
using m Ue < 5eV [^] and m„ M < 170keV p6| , respectively. Due to the much tighter 



bounds these constraints exclude all interesting regions of parameter space with such 
(Hi : /i 2 : 1^3) ratios. 



The negative sign for the neutrino mass in Eqns. (|3.3|) and (3.4) can be removed by 



redefining the fields in Eqn. ( |3.1| ) (see [^7] for more details). 



If more than one of the /ij's is non-zero, then the massive neutrino will be an admixture 
of the three neutrino basis states plus (especially for large //j's) the two gaugino and two 
higgsino states. Using the matrix Uo to diagonalize M.^, the eigenvalues are 

U$M„U = diag{M nl , M n2 , M n3 , M ni , 0, 0, m U5 } , (3.6) 



3 This factor was also obtained in [23], but additional assumptions led to compensating factors 
that canceled the suppression. 
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where the massive neutrino is defined to be the lightest of the massive states and the four 
heavier states are termed neutralinosfj Further, note that there is no reason to expect 
alignment between the neutrino eigenstates and the charged lepton mass eigenstates (e, /x, 
and r). Neutrino mass bounds in this more general case are more complicated, and better 
constraints are in fact obtainable from analysis of charged current processes (to be discussed 
in detail in the next section) than from direct mass bounds on v T (or z/ M ) [|30| . Based on 



such an analysis, Bottini et al. |3T| gave a general neutrino mass bound of 149 MeV for a 
massive neutrino that was an admixture of v e , and v T . The present case differs from 
theirs since in the RPV framework the massive neutrino can also have gaugino and higgsino 
contributions. Nonetheless, the 149 MeV mass bound is in fact applicable as will be justified 
in the following section treating charged current interactions. This bound is shown in Fig. 
6a as the upper horizontal line, but not shown for the tan/5 = 45 case depicted in Fig. 6b 
since the bound is never reached. 

Cosmological neutrino mass bounds also exist. These are usually far more stringent 
than the neutrino mass bounds discussed thus far, and upper limits of m u <35eV have been 



given [32]. However, additional assumptions about cosmology enter when determining these 



values, which are also sensitive to the decay modes of the massive neutrino, which is expected 



to be unstable. Due to these loopholes, a MeV neutrino is not cosmologically taboo |33 . 

Upper bounds on /i 5 (obtained from exact numerical calculations) throughout the /x -M 2 
MSSM parameter space are shown for tan/3 = 2, 45 in Fig. 7. Figs. 7ab show results assuming 
a mass bound of 18.2 MeV — i.e., assuming (/i x : /j 2 : /i 3 ) = (0:0:1) (these update the 
plots of || which assumed a mass bound of 24 MeV for v T ); and Figs. 7cd show results with 
the more general bound of 149 MeV. Bounds on /i 5 weaken as either M 2 and/or |/x | increase. 
As noted in ||, for high tan/3, /i 5 values in the hundreds of GeV are permitted by the v T 
mass bound. 



B. Invisible Z°- width 

The couplings of the seven Majorana mass eigenstates to the Z° boson are given by 



4 Going beyond tree level will give small masses to the two zero mass neutrino eigenstates found 
here. Trilinear RPV terms will contribute to these corrections. Attempts to fit sub-eV mass 



neutrinos such as suggested by results of the Super-Kamiokande experiment [28| into the RPV 



framework will require knowledge of these corrections. This is beyond the scope of the present 



work. See [29] for more details. 
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£ x o x o z o = - — — xhA&mCcd ~ l5$teC cd )xd Z>M , (3-7) 

4r COS U yj 

where C c d = C\ d = {UqT z uI) C( i and T z = diag(0, 0, —1, 1, 1, 1, 1) (again adopting the 
notation of [|ll|]). The partial Z° decay width into a pair of neutral fermions (with M z > 
M Xc + M Xd ) is then given by 



r ( z°^M) = - a2A 



2 I _Lv_ 



12cos 2 #„ 



c, 



erf 



(3.8) 



In the SM, the invisible Z° width, rf™ , is given by the sum of the partial decay widths of 
the Z° into the massless neutrinos. From this width, the number of SM neutrino flavors has 
been measured to be 3.09 ±0.13 In the i?-parity conserving MSSM, the decay of the Z° 
into a pair of the stable lightest neutralinos (taken as the LSP's) should also be included^, 
meaning that F^ SM > rf^. When i?-parity violation is permitted, the situation becomes 
more complicated. T(Z° — > u 5 u B ) is suppressed by kinematic factors for the now massive 
neutrino. The Z° coupling may also change due to gaugino and higgsino contributions to 
u B . Further, whether or not the partial decay width to the massive neutrino or to any of 
the more massive neutralinos should be included in will depend on how these particles 

Ainu 1 r 

decay.Q As a result, need not be larger than F§ M . An example of this is shown in 

Fig. 8, where is plotted for the generic parameter space point M 2 = fi = 200 GeV 

and tan/? = 2 assuming contributions from the two massless neutrinos and u 5 (dotted curve) 
and from the two massless neutrinos, v B , and the lightest of the neutralinos (dashed curve). 
Since the LEP measurement of Tz^ given in Table I is slightly below the SM value, the RPV 
framework can actually give better agreement. The solid horizontal lines in Fig. 8 give the 
3<T bounds for the LEP measurement. These do not pose a strong constraint on /i 5 , allowing 
values up to ~ 100 GeV in this compared to the neutrino mass constraint which from 



5 Other contributions may also have to be included — such as Z° decays to the second lightest 
neutralino which in turn decays to the lightest neutralino and a pair of neutrinos; or Z° decays to 
sparticles which are close in mass to the LSP and thus produce decay products too soft to detect 
when they in turn decay to the LSP. Lhese effects will always increase . 

6 Light scalar states might also contribute to the L^j 7 ', either directly or as virtual propagators. 
Lhis further complicates the situation due to the larger number of free parameters in the RPV 
scalar sector. We make the reasonable assumption here that such scalar states are too heavy to 
make meaningful contributions. 
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Fig. 7c demands yU 5 <21.5 GeV for m„ 5 < 149 MeV. A thorough study of the decay modes of 
the massive neutrino as well as the other particle states in the model is necessary to make 
more definitive predictions. Such input will also be required if we wish to constrain this 
minimal RPV model using the full Z° width and the searches for anomalous Z° decays, 
such as to neutralinos (Z° — > XiXjiX^ u 'i3 !)• 



IV. CHARGED CURRENT INTERACTIONS 



Up to this point, color-singlet fermion interactions with on- and off-shell Z° bosons have 
been analyzed. Important constraints can also be obtained from charged-current processes 
such as the decays of pions, leptons, and heavy neutrinos. The relevant W ± interactions 



may be written (again following notation from [ll]) as 



Wx'x — 92 TT/-U- 



\/2 



Xal» 



+ h.c. 



(4.1) 



in four-component mass eigenbasis notation. The B L ' R matrices giving the effective coupling 
strength among the mass eigenstates can be obtained from the diagonalizing matrices of the 
charged and neutral fermions: 



where 



j.L 



B L = UlT L U 



Wo ^ 



and 



B R = UlT R Un 



(4.2) 



L 



and 



T 



J 



-y/2 

1 4v 



(4.3) 



/ 



Of particular interest are pion and lepton decays. Here discussion is limited to tree-level 
decays mediated by a virtual W boson. (Z° exchange in lepton decays is negligible for the 
cases we consider below, the corresponding amplitude being proportional to the product of 
two /Vs.) An exhaustive discussion would in principle also require consideration of pos- 
sible virtual scalar intermediate states. These (together with radiative corrections) would 
re-introduce trilinear RPV terms; however, as already noted, since supersymmetric scalar 
particles must be considerably more massive than the gauge bosons, their contributions can 
be expected to be very small. Partial decay widths for these processes are then given by 



and 







7r£ 



8tt 

1927T 3 



c=l 



c,d=l 



(4.4) 
(4.5) 



where 



B ac = \ B ac 



L |2 



+ \B, 



R\2 



(4.6) 
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Standard parameters are the Fermi constant G (see [31] for minor subtleties), the pion decay 
constant /„., and V^f M = the w<i-component of the CKM quark-sector mixing matrix. TZ n i 
and TZp are leading radiative corrections to the processes. The former depends on the pion 
and charged lepton masses; the latter depends only on the mass of the decaying lepton (See 



25| , |3~4|| for more details.). Finally, the functions P^ 1 and include the entire phase space 



factors for the decays as well as parts of the matrix elements — including all dependence 
on neutrino masses. Explicit formulae from these functions are given in [31 1 ; for use below, 
Pf is: 

Pf = 9{m„ -m t - m Vc ) + 5^ - - OW 1 - %*> O ! ( 4 - 7 ) 

OtlT 1 "1/ C 7T 



As noted previously, Ref. [[H]] analyzed three-neutrino mixing. The present case differs 
since mixing of neutrinos with neutralinos (and charged leptons with charginos) is also 
possible. However, the above formulae remain valid since the new states would be too heavy 
to contribute to these GeV or sub-GeV decays. 

The experimental inputs are the following decay rates: 

r(vr -> ev e ) , r(vr -> (jlv^) , 
T(fjL-^ ) , r(r -> ez/ e z/ T ) , T(r -> nu^v T ) . 

Here, by v e ,v^v T we mean the states produced alongside the e,/z, r leptons in the corre- 
sponding decays. To eliminate the uncertainty in some of the common factors and thus 
better isolate the effects of non-SM leptonic masses and mixings, it is preferable to work 
with ratios 0,0. For pion decays, we will use 



For the lepton decays, reduced decay widths are defined by 

r " - "> = E ^ c <^ , (4-9) 

U m l' c,d=l 

which conveniently allows us to define the ratios: 

p/xe 7? v R 2 R 2 P^ e T re V R 2 R 2 P re 

fl£ = = '^^d^Peu^cd ^ ^ r e = 1 = ^d^vge^ (4JQ) 

re pre t, ^ R 2 D2 pre r M po R2 p T M v y 

1 Kt l^c,d^TU c n ev d r cd 1 l^c,d C TU c ^^u d r cd 

Barring some miraculous alignment, the physical neutrino mass eigenstates will not be 
the partners of the charged leptons — the u/s above. Nor will they be basis states in the 
single- VEV parametrization — the z/j's of Eqn. fl3.1|). The physical eigenstate which acquires 
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a mass at tree- level is denoted by u s ; we will denote the other two massless degenerate 
states by u zl and u z2 . Note that these two light eigenstates are not uniquely defined at tree- 
level; any linear combination between them is another massless neutrino eigenstate. This 
indeterminacy will not affect physical results; moreover, the degeneracy will be lifted by 
radiative corrections. 

In the rest of this section, we will discuss implications of experimental constraints on 
the ratios (|4.8| ) and A4.10|) . Exact numerical values for rotation matrices Ul,Ur, and U Q 



will be used to compute the charged current couplings B L and B R . However, in order to be 
able to understand the qualitative features of our numerical results, we find it worthwhile to 
derive and use approximate analytical expressions for these couplings. To this end, we will 
use the following approximation. We split the diagonalizing rotation matrices into separate 
"quasi-MSSM" and "quasi-SM" blocks: 

U L = diag-f/^, J 3x3 } , U R = diag{i? fl , J 3x3 } , and U = diag{[/ 4 , R 5 } , 

where U4 diagonalizes the first 4x4 block of M. M and -R5 is again the 3x3 CKM-like 
neutrino mixing matrix. Note that, while this approximation is valid in the small /x, limit, 
it is not the same as the approximation used in section II. In that case, all first order terms 
in the Hi were kept; here, we keep only those terms of the form fii/ '^5. This will be the 
meaning of the "small fi approximation" in this section. 

The SM components of the z/ 5 eigenstate are then given by 

Rf = - , (4.11) 

Ms 

and, by unitarity, the quasi-SM components of the other two eigenstates obey the sum rule 

\I$\ 2 + \I%\ 2 = l-(**\ . (4.12) 



Using Eqns. Q4.11| ) and Q4.12| ) together with the assumption of no right-handed neutrino 



fields, the quasi-SM sector charged current couplings become 

B R = B L = J{ mn > B 2 = \R mn \ 2 (4 13) 



A. Pion decays 



Inserting the small- /ij approximate expressions from Eqns. ( 4.11| ), ( |4.12j ) and ( 4.13 ) into 

contributions involving the more massive quasi-MSSM 



Eqn. ( }4.8|) , and disregarding -B| mX 
neutral states (as justified earlier), yields (using PJ[ e = P£ e for the massless states): 
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(4.15) 



is the SM prediction (with = 0.96103), and P rf = P* e /P z f - 1 . 

The behavior of the i?^ ratio as a function of /is can be understood by analyzing the 
dependence of the kinematic functions P B on the neutrino mass. Indeed, for a given set of fii 
ratios, the only dependence on /is in comes via m U5 in Pf. Small chan e;es in m vc . affect 
P B e and P 5 nfJ- differently. P 5 ne increases quite rapidly with increasing m U5 , as a consequence 
of matrix element dependence on the masses of decay products in pion decay [for m M5 < 80 
MeV, P B e = (m U5 /m 7T )' 2 is a very good approximation]. In contrast, P^ decreases, albeit 
more slowly, because of the phase space factor. As a result, their ratio R™ increases with 



m 



and this effect is further enhanced the larger '/x 5 is with respect to ^2/^5 ■ 



Looking in more detail at various mass ranges, we find that for m^ 5 <20 MeV, the kine- 
matic function P^ is nearly the same as the constant PJl^ (V nfl ~ 0). In this case, the 
experimental constraints on the ratio R™ can be put in a simple form in terms of the RPV 
parameters. Using the 3cx bounds from Table I, we obtain: 



for 



m UB <20 MeV 



/'1 



< 



/'5 



137 m h 



(4.16) 



- m, 



Beyond the threshold value of m U5 
the massive neutrino can no longer proceed; P^ p 
constant. Then we have the following constraint on the ^2/^5 ratio 

for 34 MeV < m U5 < 139 MeV 



33.91 MeV, the 7r decay into a muon and 
= and the denominator of (|4. 14|) becomes 



1*2 



< 



/'5 



137 



(/i! = 0). 



(4.17) 



Upper limits on the /X1//X5 ratio in this region are given in Fig. 9. 

For values of m U5 of order 100 MeV, the increase in P B e due to the amplitude contribution 
are offset by the decrease due to the phase space factor, and R™ begins to decrease. At 
the threshold m U5 = — m e = 139.057 MeV, the ir cannot decay into ev^ either. Above 
this threshold, the R™ ratio is constant (in the approximation used here), and the following 
constraints can be derived: 



m U5 > 139 MeV 



for ix 2 > Mi 
for > /i 2 



i4 f4 ^ 1 

A? A? < 137 
,,2 



— 7" — — T - 



1 

82 



(upper experimental bound), 
(lower experimental bound). 



(4.18) 
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It must be mentioned that large neutrino masses of (9(100 MeV) may be beyond the range 
of validity of the small-/ij approximation unless tan/3 is "small". For tan/3 = 2 and for the 
MSSM parameters M 2 = 250 GeV and fj, = 100 GeV, Eqn. (0) gives m Uh = 149 MeV when 
/i 5 = 13.9 GeV, in rough agreement with the see-saw prediction of fi 5 = 11.4 GeV from Eqn. 
( p.3| ), and thus arguably still within the small fi 5 domain. In addition, radiative corrections 
to this tree- level calculation may alter the numerical values given above. 

Using exact numerical calculations, Fig. 10 shows the dependence of R™ on /i 5 for some 
illustrative ratios among the /Vs and for the generic MSSM parameter point M 2 = fi = 
200 GeV and tan/3 = 2, 45. The solid horizonal line denotes a 3cr positive deviation from the 
experimentally determined central value. The qualitative attributes of the curves agree well 
with the predictions made above on the basis of the analytic expression flOp . Note that 
for large tan/3, high values of m U5 are unattainable no matter how large /i 5 is. In this case, 
values of \i h in the hundreds of GeV are allowed. 



B. Decays of charged leptons 



Eqn. ( |4.5|) shows that H^*-mediated charged lepton decays are proportional to two B\ aVc 
factors. This along with the more complicated P^/ function (which gives the dependence on 
the neutrino mass) makes analytic expressions unwieldy, even in the small-/Zj approximation. 
Hence only results from the exact numerical calculations will be discussed. 

R^l of Eqns. ( [4.1 0|) provides the most interesting constraint, tightening the bound on ^ 



in some regions where is less effective. A similar conclusion was reached in [31]. For 



tan/3 = 2 and /i x set to zero, exact numerical results at the generic MSSM point M 2 = /x = 
200 GeV show R^l to be more restrictive than R^ for ^ < and \ < < 1 (echoing 
Fig. 6 of pi]). Fig. 11a depicts the actual behavior of R^ e e vs. /x 5 at this M 2 -/i point for 
tan/? = 2 and assorted /ij ratios. The solid horizontal lines represent 3a deviations above 
and below the experimentally-determined central value (which is consistent with the SM 
prediction of 1). For /i 2 //i 3 closer to one (g^^ < 1) and smaller values of /i 5 (yU 5 <10 GeV), 
R^l may dip below the experimentally-allowed band; however, for any n % j \i z < 1, R%% 
eventually goes above the acceptable band as \i h is increased. The upper bound on /i 5 only 
runs from 19.5 GeV down to 17.5 GeV as ^ is changed from i to (values below are 
indistinguishable from 0). However, the neutrino decay constraints from the WA66 and 
CHARM experiments described in the next subsection may surpass the constraint from R^ e e 
for these somewhat larger /x 5 values. 

For high tan/3, the behavior of R^ vs. /i 5 is decidedly different, as seen in Fig. lib for 
tan/3 = 45. Here the value of R^ e e always drops below the experimentally- allowed region as 
/i 5 is increased irrespective of the /i 2 : fi 3 ratio. The behavior is qualitatively reminiscent 
of that for the ^ = 1 curve for tan/3 = 2 even if now ^ is set to zero. Quantitatively 
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though, the upper limit placed on u 5 becomes quite large as — drops to zero. This bound 
is nevertheless still more stringent than that from R 1 ^ if fj, 1 ~ 0, again irrespective of 

The experimentally-derived values for T Me , Y re , and T rAi can also be applied individually 
without taking ratios. The 3a bounds on these quantities are also given in Table I. However, 
these restrictions were always found to be weaker than the constraints from the ratios (in 
contrast to Ref . |3T| ) . 



C. Decays of a massive neutrino 

Assuming the decay of the massive neutrino is mediated by a virtual H^-boson, the 
expected decay modes are v\ — > W^P; W*^ 1 — > ^v\ t , (qq f n^. 

A crucial experiment restricting this process was performed by the CERN WA66 Col- 
laboration |^5| using the BEBC bubble chamber placed in a neutrino beam resulting from 



dumping protons on a high density target. The neutrino beam is mainly composed of 
'prompt' neutrinos from charmed meson decays which can include massive neutrinos up to 
~1.8GeV. WA66 is sensitive to massive neutrino decays into electrons, muons, and pions. 
From the absence of any excess of such events, limits can be placed on mixing of the massive 
neutrino state with either the u e or weak-flavor eigenstate. In the smail-/Zj approximation, 
vi = Vi and this mixing is simply |_R™ n | 2 = (jj^J . Therefore, the WA66 results can be used 
to restrict these ratios if the small-/Zj approximation is valid; for m Vh ~ (9(100 MeV), this 
would demand that tan/5 is small. 

The WA66 results can be summarized as follows (results from CHARM |36| are similar 
in the parameter regions of interest here). For neutrino mass values of order 100 MeV and 
above, the mixing parameter \x\j \x§ has to be smaller than 10~ 3 (< 10~ 4 , for m„ 5 ~ 149 
MeV). The constraint on \x§ is a little bit weaker: at m UB = 149 MeV, < 10~ 3 . 

These limits as a function of m V5 are presented in Fig. 9. Note that for m U5 less than 0(80- 
100) MeV, the constraint on nif ^5 coming from 7r decay is stronger than that from WA66, 
while for m U5 > (9(80-100) MeV, the WA66 constraint is stronger. Thus, these two inputs 
play complementary roles in setting limits on neutrino flavor-state mixings. 

The WA66 neutrino decay limits are also used in establishing for the present RPV sce- 
nario the 149 MeV absolute upper mass bound on a neutrino which is not a pure z/ r |] As 



7 This result from Ref. [31| is based on a charged current analysis of R^, R T T e ^ Rlj%, T^ 6 , T Te , 



and T Tfl . As mentioned earlier, Ref. [31] only includes 3- flavor mixing, not possible mixing with 
gauginos and higgsinos. It should also be noted that this value is given in |3l| as a la bound 
whereas more conservative 3a bounds are used throughout the present work. 
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seen above, a neutrino with a mass of 149 MeV is about 99.9% u 3 which is now extremely 
well aligned with the r lepton. In this case the tighter LEP bound of m Vr < 18.2 MeV is 
applicable; therefore, the upper mass limit on a mixed massive neutrino state is ~149MeV 
as adopted here. 

There are other experiments on neutrino interactions which constrain the mixing among 
the three generations (for example, the CHARM II experiment [57] on — e scattering, 



and LAMPF |38fl on v e — e scattering). Some of the resulting constraints are stronger than 
the WA66 results. However, the interpretation of these data in our current framework 
(constraints on the RPV parameters /ij) is straightforward only in the approximation used 
in this section; otherwise, mixing between the neutrino and neutralino contributions should 
also be taken into account. For this reason, we have discussed only the WA66 results (for 
purposes of illustration), leaving a more detailed analysis to a future paper. 



V. NEUTRINOLESS DOUBLE-BETA DECAY 

Neutrinoless double beta decay (Of/3/3) places very stringent bounds on trilinear RPV 
parameters, as demonstrated in Ref. fl39"f . The bilinear /i« RPV couplings can also mediate 
Ov/3/3] therefore, this process warrants serious consideration here. First, however, features 
unique to this process must be carefully noted. As with pion decay, this process involves 
(valence) quarks rather than leptons in the initial state. However, unlike pion decay, for 
QvfiP the decay products are not all colorless. Thus transitions can be mediated by W- 
bosons and a Majorana neutrino with mixing to v e or by scalars (squarks and/or sleptons) 
and sfermions (gluinos or neutralinos) f40| , |39"|j . Tree- level diagrams involving the former pair 



can derive the neutrino couplings from the RPV bilinear /ij's, while those involving the 
latter pair will be proportional to trilinear RPV couplings. For leptonic decays, possible 
tree-level trilinear RPV A coupling dependence due to scalar intermediates was neglected on 
the grounds that this would be kinematically suppressed since sleptons and Higgs bosons 
are more massive than the gauge bosons (and the trilinear couplings are expected to also be 
small). In contrast, with Qvf3f3 there are strongly-interacting sparticle intermediates which 
could perhaps off-set the simple kinematic suppression; thus, trilinear RPV couplings might 
play a significant role at tree-level for this process. It should be stressed that this is not 
the case for the previous sections — there the trilinear RPV couplings have not been set to 
zero by hand, they simply do not yield any significant contributions to those processes at 
tree-level within the single- VEV parametrization augmented by the reasonable assumption 
of sufficiently heavy scalars. For OvP/3, in order to concentrate on the bilinear RPV /Vs, 
the A' couplings can be assumed to be negligible and/or the squark and gluino masses can 
be assumed to be very large to kill these contributions. Alternatively, to set conservative 
bounds on the //j's it is sufficient to assume that there is no destructive interference between 
the two types of diagrams. A more thorough analysis including the scalar intermediates is 
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underway. 

With these caveats, the effective constraint from 0uj3(3 becomes EI 



™**\ B evr\ < 0.46 eV for m U5 < lOMeV . (5.1) 



Applying the small- /ij approximation, this translates into 

/046 




x 10~ 3 (with m us in MeV). (5.2) 



An alternative approximation is to consider only the leading mixing effect between the sole 
vev-bearing Y = | superfield basis state and the other Y = | superfield basis states. This 
yields the approximate expression [cf. Eqn. ( [4.131) 1 



and thus \& „ , | 2 + \B\ v f ~ 1 - ^ 2< 4+/o) 



(5.3) 



This approximation should hold even for large values of /z 5 . Using this approximation along 
with Eqn. (j3.4|) for m^, the 0^/3/3 constraint is 



^ < K 2.15 x 10-7(1 (5.4) 

//„ f cos/3 \{xg% + g() v 

(with mass parameters in GeV). For the low neutrino mass region in which it is effective, 
the 0vf3(3 constraint on // a (or ^) is much stronger than that from pion decay, as is seen in 
Fig. 9. At the generic MSSM parameter point (M 2 = //„ = 200 GeV), Eqns. (^) and ( ggfl 
require ^< ^ for MeV scale neutrino masses and : /i of around 1 : 150, 000 (1 : 15, 000) 
for tan/3 = 2 (45) — setting ~ a U° ws Me values of only tenths of a GeV for tan/3 = 2 
and a couple GeV for tan/3 = 45. 

The limitation from Oz//3/3 is only applicable if m U5 < 10 MeV. Thus the Ov/3/3 injunction 
is turned off if fi 5 is large enough to push m VB above this threshold. This means that at any 
given point in MSSM parameter space, and for fi ± ^ 0, there will be at least two allowed 
ranges for /x 5 : /i 5 < fj,P< max and < /i 5 < /j™ ax , where ^ max is the upper bound from 
Oz^/3/3, /ii 10 ' ) is the /i 5 value at which m U5 reaches 10 MeV, and /i™ ax is the cut-off value due 
to the strongest constraint aside from Oz//3/3 (typically this is from R™ for low to moderate 
tan/3 values and from Rl^l for high tan/3). This of course assumes ^ max < /4 10 - 1 < jj^ ax as 
is almost universally true. 



VI. OVERALL COMBINED CONSTRAINTS 

Here we pull together all the constraints addressed individually in the preceding sections. 
These are combined numerically in a comprehensive program to yield a maximum- allowed 
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/x 5 for any given point in MSSM parameter space for a specified fi ± : /x 2 : /x 3 set. Table I lists 
all the experimental constraints applied. Note that constraints resulting from the WA66 and 
CHARM neutrino decay experiments as well as those from neutrinoless double beta decay 
experiments are not implemented except to rule out potentially admissible large /X5 regions 
beyond the first cutoff point for individual constraints (a thorough implementation awaits a 
more complete analysis including scalar intermediates). 

A few details associated with the numerical studies deserve mention. First, (electroweak) 
gaugino unification is assumed; i.e., M ± = xM 2 with x = ^ta,i?9 w . Second, running parame- 
ters are evaluated at the scale appropriate to each particular process. Thus for instance Z° 
decays use sin 2 9 W (M Z ) while r decays have sin 2 ^„,(m T ). Specific numerical inputs include 
(see 0): sm 2 6 w (M z ) = 0.2315, M z = 91.1867 GeV, and M w = 80.4 GeV. 

Fig. 12 scans the M 2 -/x plane and presents contours of maximum- allowed /x 5 values in 
GeV. Figs. 12a and 12b have tan/3 = 2 while 12c and 12d have tan/3 = 45; /x x : /x 2 : 
/x 3 = : 1 : 1 in 12a and 12c whereas fi 1 : /x 2 : /i 3 = : 1 : 10 in 12b and 12d. With 
tan/3 = 2, the R^% constraint dominates for the 0:1:1 ratio set; but R 7 ^ is stronger for the 
: 1 : 10 ratio set. The charged current constraints are slightly more restrictive than those 
from the 149 MeV neutrino mass bound shown in Fig. 7c — the difference is greater when 
/Xi : /x 2 : /x 3 = : 1 : 10, but the order of magnitude is still the same. If /x 2 is set to zero 
(fjL-i : /x 2 '■ fa — : : 1), then the 18.2 MeV v T mass bound of Fig. 7a dominates. 

Another distinction between Figs. 7 and 12 is in the central region where M 2 or |/x | 
is small. This region is ruled out in Fig. 12 by limits from processes involving charginos 
and neutralinos. Such restrictions, embodied in the Table I bounds on the chargino mass, 
on anomalous "visible" Z° decay modes (Z° — > x ± £ =F , x Q c Xdi Xd v \ d 7^ 1)> an d on Z° full 
and invisible decay widths, are in general rather conservative, reflecting present intangibles 
concerning sparticle decays — a more complete study of these is in progress. As discussed in 
Section II. E. and seen in Fig. 1, /Xj values sufficient to significantly affect the lighter chargino 
mass are in the high-/Xj realm and thus only depend on /x 5 . Results for the combined 
constraints concur: bounds in the small M 2 , small |/i | region are basically independent 
of the ratios among the /Vs. One additional point worth noting is that driving up /x 5 to 
extremely high values will not be able to push up a very small MSSM chargino mass. 

Regarding tan/3 = 45: the general pattern of the contour lines is similar to that for 
the tan/3 = 2 plots; however, admissible /x 5 values are much larger. Upper bounds on 
/x 5 are relaxed by a factor of ~ tan/3. The limits in Figs. 12c and 12d are much stronger 
than those in Fig. 7d which just uses the 149 MeV neutrino mass bound. Both of these 
changes are due to the dominant charged current constraints (Rtf% is strongest except for 
/Xj : /i 2 : /i 3 = : 1 : 10 and M 2 <160GeV in which case R%£ dominates). As /x 2 — > 0, we 
again return to the 18.2 MeV m VT limits of Fig. 7b (with the central region again excluded 
by chargino and neutralino constraints). 
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VII. CONCLUSIONS AND OUTLOOK 



The degree of generality possible within the single- VEV parametrization deserves special 
notice. With only a trio of inputs (the bilinear /ij couplings in the superpotential) beyond 
those of the MSSMQ, the tree-level mass matrices of all color-singlet fermions are completely 
determined, and a broad range of leptonic phenomenology can be analyzed with a reasonable 
level of sophistication. Furthermore, the trio can often be collapsed into a single input, // 5 , 
which carries the full weight of i?-parity violation. This is to be contrasted with analyses 
which either contain a plethora of RPV parameters which preclude a meaningful coverage of 
the parameter space within the model or arbitrarily pick one or two RPV paramaters to be 
non-zero while the others are all set to zero by hand. As a result, the true freedom within 
an RPV model will be masked or muddied by a less optimal choice of flavor basis (or by no 
choice at all!). 

We have presented tree-level coverage for the variety of electroweak signals summarized 
in Table I. In the single- VEV parametrization analysis presented here, the trilinear RPV 
couplings have not been set to zero. Rather, they simply do not contribute (at tree level 
with gauge boson propagators) to the impressive spectrum of leptonic processes studied. 
In principle, intermediate scalars (sleptons, Higgs bosons) can usher back in the tree-level 
trilinear dependence; however, such contributions should be suppressed by the larger masses 
of the scalars relative to the SM intermediate vector bosons, and are nonnegligible only in 
special regions of the parameter space. A more in depth look at the scalar sector (along 
with the related issue of electroweak sparticle decays) is now in progress. An exception to 
ignoring the scalar intermediates may be necessary with neutrinoless double beta decay due 
to the presence of strongly-interacting sparticle intermediates. This question is also under 
study. 

Also beyond the present study but slated for future work are loop effects, necessary 
to describe interesting processes such as /i — > ej. This process may place an additional 
significant restriction on the /Vs. Other loop processes will be important in constraining the 
trilinear RPV couplings. Loops will also lift the degeneracy of the two tree-level massless 
neutrinos. This extra degree of precision is certainly needed to study the very low-mass 
neutrinos preferred by several neutrino oscillation experiments. Some models to describe 
such experiments also suggest that one or more extra light (sterile) neutrinos be added. 
This could strongly affect the analysis of the charged current constraints presented here, 



Of course the MSSM already has a fair number of input parameters — M 2 , fi , and tan/? all 
enter into the mass matrices. This is unchanged and precludes an exhaustive scan of the full 
beyond-the-SM parameter space. 
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but not the Z°-mediated neutral current constraints for the charged leptons. Thus, while 
the charged current constraints are more restrictive, those for neutral current processes are 
more robust against such possible model extensions. 

The single tree-level neutrino mass, m UB , depends on the three fa only through fa.. If 
MeV-scale neutrino masses are allowed, then fa may well be large enough to yield signal rates 
near, at, or above present experimental bounds for the numerous other processes described 
herein. Such a neutrino mass is not ruled out by direct or indirect machine (terrestrial) 
mass bounds. Cosmological constraints favoring light neutrinos may not be applicable; this 
depends on the decay properties of the massive neutrino and requires more study. An MeV- 
scale neutrino has been found to be consistent with cosmology in at least one study (which 
did include fields beyond those in the MSSM) 

The results of this analysis reinforce those of : fa. values of the same order of magnitude, 
or even much larger, than M 2 and fa (the standard MSSM inputs) are allowed by the 
experimental bounds. This is particularly true for high tan/3 (>45), where RPV signals are 
strongly suppressed. In this case even much tighter bounds on the neutrino mass do not 
preclude large fa values. This tan/? dependence is almost universal among the processes 
studied. This plus the fact that both the neutrino mass and numerous (though by no means 
all) approximate expressions for other bounds depend only on fa lead to the interesting 
question of how much variation is possible among these other processes if fa is fixedQ. To 
look at this various possible fa : fa : fa combinations were studied using the exact (at 
tree-level) numerical expressions. Q A rough hierarchy is seen in the constraints on the 
individual fa's: fa is strongly restricted, fa is less restricted, and fa is still less restricted. 
This suggests ratio sets of the general form : 1 : x (x > 1). With only two free inputs, this 
further suggests that, along with a dimensionless ratio, fa is the preferred indicator of RPV 
effects in leptonic phenomenology. 
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9 If both fa and neutrino mass are fixed, then restrictions are placed upon the MSSM input 
parameters. How much tan/? dependence and variation is then possible in the remaining processes 
is currently being investigated. 

10 Naturally, studies that truncate the number of generations cannot preform such analyses and 
miss very significant and interesting effects. 
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Table Captions 



Table I: 

Summary of phenomenological constraints incorporated in the overall parameter con- 
straint plots. The invisible Z°-width is assumed to include decays to end states composed 
of neutrinos and the lightest neutralino. Bounds imposed on Z° decay constraints involving 
charginos and neutralinos are quite conservative, representing uncertainty in signal detec- 
tion. 



Table II: 

Experimental left-right asymmetry results: Ai obtained from A e FB (using A l FB = \A e Ai) 
measured at LEP and SLD, V T measured at LEP, and a direct measurement from SLD. 
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Figure Captions 



Figure 1: 

Branching ratio for — > e~e + e~ as a function of \x h (in GeV), with M 2 = fi = 200 GeV 
and tan/3 = 2 (left) or tan/? = 45 (right) for ratios of \i x : /i 2 : /i 3 as marked. The solid 
horizontal line is the experimental bound. If fi 3 — here set to zero — is varied, it affects 
the curves only via /i 5 , and hence stretches out the horizontal scale. 



Figure 2: 

Branching ratio for t~ — > e~ [i + as a function of fi 5 (in GeV), with M 2 = fi = 
200 GeV and tan/3 = 2, for ratios of [i± : /i 2 : /i 3 as marked. The solid horizontal line is 
the experimental bound. Here, /x 2 is set to zero; again, varying this stretches the horizontal 
scale. 



Figure 3: 

Leptonic partial decay widths of the Z° (in MeV) as a function of /i 5 (in GeV), with 
Hx : /i 2 : fi 3 = 1 : 1 : 1 (so that \ii = -/j/-^); M 2 = fi = 200 GeV, and tan/3 = 2. Horizontal 
lines at the left edge of the plot are the ±3<r experimental bounds |yj for the corresponding 
quantities; they are independent of yU 5 and are truncated at right for clarity. Explicit values 
for the experimental bounds on these partial decay widths are found in |12| . 



Figure 4: 

Leptonic L-R asymmetry: deviations from SM predictions as functions of fj, s (in GeV), 
with jj, x : \i % : fi 3 — 1 : 1 : 1 and M 2 = fi = 200 GeV for tan/3 = 2 (top) and tan/3 = 45 
(bottom). The solid line approximates the experimental 3a upper bounds for each of the 3 
asymmetries. 



Figure 5: 

Chargino masses. Contours show minimum values of /i 5 (in GeV) necessary to push the 
lighter chargino mass above 90 GeV for tan/3 = 2 (top) and tan/3 = 45 (bottom). The region 
above or outside a given contour has M cl > 90 GeV for /i 5 at or above the designated value. 



Figure 6: 

Neutrino mass m U5 as a function of fi 5 (in GeV), with M 2 = /j, — 200 GeV and a) 
tan/3 = 2, b) tan/3 = 45. The lower horizontal line is the 18.2 MeV machine bound for a 
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pure u T ; the upper horizontal line is the 149 MeV bound for a generic u 5 . (M 1 = xM 2 , with 
x = ^tai?9 w .) Insets: Low /i 5 portion of curve on log-log scale. 



Figure 7: 

Maximum allowed values of /x 5 (in GeV) consistent with neutrino mass bounds: machine 
bound m UT < 18.2 MeV (applicable for fi 1 : /i 2 : /x 3 = : : 1) for a) tan/5 = 2 and b) 
tan/? = 45; the absolute bound m„ 5 < 149 MeV for c) tan/3 = 2 and d) tan/3 = 45. The 
region below or inside of a given contour is excluded for /x 5 's above the indicated value. 



Figure 8: 

The invisible Z°-width (in MeV) as a function of \i h (in GeV) , with M 2 = // = 200 GeV 
and tan/3 = 2. The solid horizontal lines are the upper and lower experimental bounds. 

is assumed to be T(Z° — > u c Ud) (dotted curve) or T(Z° — > v c Vd,v c X®,XiXi) (dashed 
curve) where Xi is the lightest neutralino (4 th lightest neutral color-singlet fermion). 



Figure 9: Constraints on and /i 2 / 7*5 ratios. Solid lines: pion decay [cf. Eqn. 

( |4.14 )]; in the upper half of the figure, the lower solid line corresponds to ^2/^5 — 1/12, 
while the upper line to ^2/^5 = 0; in the lower half of the figure, the solid line corresponds to 
/^i/^s = 0. Dotted lines: the WA66 experiment [cf. R35H 1. Dashed line: neutrinoless double 
beta decay (cf. Eqn. ( |5.2j )). 



Figure 10: 

K* e as a function of /x 5 (in GeV) for various m ratios, with M 2 — fi — 200 GeV and a) 
tan/3 = 2, b) tan/3 = 45. In a) the /ij ratios are 1:1:1 (long-dashed line); 1 : 10 : 110 
(dot-dashed line; appears twice in the figure); 0:1:1 (open-spaced dotted line); : 1 : 11 
(closely-spaced dotted line); 1 : : 110 (short-dashed line). 



Figure 11: 

vs. fj, 5 for assorted /ij ratios and a) tan/3 = 2, b) tan/3 = 45. Horizontal lines denote 
3cr deviations from the measured central value as given in |3l] . 



Figure 12: 

Maximum allowed values of /i 5 (in GeV) consistent with all the constraints listed in Table 
I: a) for tan/3 = 2 and fii ratios : 1 : 1; b) for tan/3 = 2 and : 1 : 10; c) for tan/3 = 45 and 
: 1 : 1; d) for tan/3 = 45 and : 1 : 10. 
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TABLES 



TABLE I. 



Quantity 



Hi combo, 
constrained 



Experimental bounds H 



Z -coupling: 



K T 

Br(Z° 
Br(Z° 
Br(Z° 

Br(fi~ 
Br(r~ 
Br(r~ 
Br(r~ 
Br(r~ 
Br(r~ 
Br(j- 
Br(Z° 
Br(Z° 
Br(Z° 

r z 

^ mi > 



(e-fj, universality) 
(e-T universality) 
(fj,-T universality) 
(e-/i L-R asymmetry) 
(r-e L-R asymmetry) 
(t-/j> L-R asymmetry) 

e±rT) 
-/iV^) ^ 
— > e e + e ) 
— > e~e + e~) 
-> /i~e + e~) 
-» /i + e~e~) 
-► e-fi+fi-) 
-+ e + fi-fi-) 

-* M~M + M~) 
- X^ T ) 
X ± X T ) 



o,,o 



(total Z°-width) 
(invisible Z° width: 

-> VcVd,V c Xl,XlXl) 



IT^-coupling: 










ez/i^) 


.r e 


(r-> 


ez^) 




(r-> 


jJLVv) 




(-7T decays) 




(r decays) 




(decays to e's) 


• "V 5 


DL |2 




mass constraints: 




• z^ 5 mass 






• x^ 1 mass 







Ms 
Ms 



Mi 
M2 



M? -M2 
M? "Ms 
M2 - Ms 
Mi — M2 
+ Rt. contrib. 
+ Rt. contrib. 
IM1M2I 

ImiMsI 

IM2M3I 

IM1M2I 

ImiMsI 

IM2M3I 
Im?m 2 msI 

ImiMsI 
IM1M2M3I 

IM2M3I 

Ms 
Ms 
Ms 
Ms 
Ms 



m^ 5 / /ij ratio 
m^ 5 / /Xj ratio 
m Us I fa ratio 
m /ML aTlr l M 

m^ 5 / fii ratio 
m„ 5 / fj>i ratio 



Ma 
Ms 
Ms 



(0.596 ± 4.37) x 10~ 3 

(0.955 ± 4.98) x 10~ 3 

(1.55 ±5.60) x 10~ 3 

(0.346 ± 2.54) x 10" 2 (from Ug) 

0.0043 ±0.104 

0.082 ±0.25 

< 1.7 x 10~ 6 

< 9.8 x 10~ 6 

< 1.2 x 10~ 5 

< 1.0 x 10" 12 

< 2.9 x 10~ 6 

< 1.7 x 10~ 6 

< 1.5 x 10~ 6 

< 1.8 x 10~ 6 

< 1.5 x 10~ 6 

< 1.9 x 10~ 6 

< 1.0 x 10" 5 

< 1.0 x 10~ 5 

< 1.0 x 10~ 5 
2.4948 ± 0.0075 GeV 
500.1 ± 5.4 MeV 



0.983 ±0.111 
0.979 ±0.111 
0.954 ±0.108 
(1.230 ±0.012) x 10~ 4 
1.0265 ±0.0222 
1.0038 ± 0.0219 

< 0.46 eV (only for m U5 < 10 MeV) 



< 18.2 MeV if v s = v T 

< 149 MeV if v h / v T 
> 70 GeV 
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TABLE II. 



An 


Method 


LEP combined \ Y1\ 


SLD [17| 




FB 


0.1461 ± .0110 


0.152 ± .012 




0.1399 ± .0073 


— 


direct 




0.1543 ± .0039 




FB 


0.1488 ± .0170 


0.102 ± .034 




FB 


0.1753 ± .0210 


0.195 ± .034 


v T 


0.1411 ± .0064 
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Figure 2 



Figure 3 
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Figure 4 
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Figure 6b 
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Figure 7 
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Figure 10a 
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Figure 11a 
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